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Geometrical constraints imposed on liiglier dimensional harmonic lattices generally lead to non- 
linear dynamical lattice models. Helical lattices obtained by such a procedure are shown to be 
described by sine- plus linear-lattice equations. The interplay between sinusoidal and quadratic 
potential terms in such models is shown to yield localized nonlinear modes identified as intrinsic 
resonant modes. 

PACS numbers: 05.45.Yv, 63.20.Pw, 03.40.Kf 



(N 



> 

in 
o 

O 



> 

X 
S3 



I. INTRODUCTION 

Recent developments in the physics of intrinsic local- 
ized modes (ILM) d have shown that the ILMs are 
not only conceptually important in theoretical and math- 
ematical physics as ubiquitous fundamental modes in 
nonlinear, discrete physical systems but also possess 
innovative potential applications in fundamental science 
and technology in addition to their experimental obser- 
vation iJl, l5i, l2j . Areas of such applications cover coupled 
Josephson junctions, photonic crystals, optical lattices in 
Bose-Einstein condensates, all-optical logic and switching 
devices, targeted breaking of chemical bonds, and so on 


Very recently, the interplay between nonlinear dynam- 
ics and geometry has attracted particular attention in 
ILM problems j^]- Historically, such an observation was 
first made when formulating model Hamiltonians for the 
dynamics of bases in DNA to take care of helical struc- 
ture, to the best of the authors' knowledge The 
latter problem gave rise to the intuitive introduction of 
the so-called sine-lattice (SL) model in which inter- 
site interactions along a given strand are taken to be 
sinusoidal rather than the conventional quadratic ones. 
The SL model later turned out to yield a novel type of 
ILM, referred to as an intrinsic rotating mode (IRM), also 
termed a roto-breathcr ^12.] . An example of transition 
from oscillation to such rotation modes in simulations of 
molecular crystals can be found in |l3l |. Another exam- 
ple of this type is given by curved or bent chains 3| 
long-range interactions on a fixed curved substrate [l5j |. 

Geometric constraints naturally appear in crystalline 
bodies consisting of relatively rigid atomic clusters. An 
important example of this class of materials are the poly- 
morphs of silica (Si02), where the structural units, i.e., 
Si04 tetrahcdra are corner- linked by oxygen atoms, and 
the energy cost of deformation of the tetrahedra is much 
greater than the cost of their mutual rotations. Atoms 



in the almost rigid clusters move as if they were sub- 
jected to a geometrical constraint. To describe the po- 
sition of a finite-size molecular cluster, one has to intro- 
duce not only translational but also rotational degrees of 
freedom. It has been demonstrated that the rotational 
degrees of freedom can be responsible for the incommen- 
surate phase in quartz 0, [l9l |2p | . negative Poisson's 
ratio of cristobalite and quartz |2lL 123. l23|.|24ll , and neg- 
ative thermal expansion of beta-quartz |25| . Similar ef- 
fects can be observed in other materials with microscopic 
rotations, such as perovskites (e.g., SrTiOs) containing 
corner-linked TiOe octahedra, the KII2PO4 (KDP) fam- 
ily of crystals with comparatively rigid PO4 tetrahedra, 
among others. 

Motivated by these investigations, we proposed in a 
previous paper a mechanical model in which a set of 
masses on a linear chain are re-arranged to slide on fixed 
rings |l6j | . Such a model was shown to contain rich non- 
linear dynamics, exhibiting various types of nonlinear 
modes ranging from kinks to IRMs. The method em- 
ployed there amounts to applying a specific geometrical 
constraint to a purely harmonic lattice, leading eventu- 
ally to equations having the form of a extended version 
of the SL equations. 

It is our purpose here to propose a method for studying 
the above-mentioned interplay between nonlinearity and 
geometry in a more general and transparent way. The 
basic point of our method is to take three-dimensional 
(3D) harmonic lattices as a starting point on which ge- 
ometrical constraints are imposed. Since linear systems 
constitute a basis for studying physics and mathematics 
in general, we expect the present method will give much 
more insight into the problem. The work is composed of 
two parts: the first presents a general scheme of geomet- 
rical constraints, while the second complements it with 
a detailed study of the ILMs in sine- plus linear- lattices 
resulting from the helical constraints that we impose on 
the original 3D harmonic lattice. 
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This paper is organized as follows: in the next sec- 
tion, we present a simple 3D harmonic lattice model and 
consider a general scheme of the geometrical constraint 
applied to it. In Section III, we introduce the helical 
constraints as an application of the general method to 
arrive at helical lattices described by sine-linear-lattice 
(SLL) equations interpolating between SL and linear lat- 
tice. Section IV is devoted to the study of some prop- 
erties of the SLL equations. Generalization of the SLL 
model is made in Section V to include the effect of an on- 
site potential, giving rise to the ILMs. The last section 
is devoted to concluding remarks. 



II. THE LINEAR LATTICE MODEL AND A 
GEOMETRIC CONSTRAINT 



We consider a three-dimensional lattice governed by 
the Lagrangian L 



K 



where T and U are the kinetic energy and the poten- 
tial energy of the system, respectively. The quantities 
Xn,yn, Zn are the dynamical variables associated with the 
nth atom of atomic mass m„. Then, the Euler-Lagrange 
equations assume the form 



<fx. 



' dt 

d^Z^ 



K{Xn+l + Xn-l - 2x„), 
2 - K{yn+l+yn-l~'2.yn), 



rrir. 



K{Zn+l + Zn-1 - 2z„). 



(2) 



Physically, Eqs. (O are equations of motion for harmonic 
lattice vibrations of a simplified version of a simple cubic 
lattice, in which a;„, ?/„, z„ represent the x— , y—, z— com- 
ponents of the displacement vector r„ of the nth atom 
from its equilibrium position. 

Suppose now that there exists a single variable s such 
that 

x„ = /(s„), yn = .9(s„), Zn = h{s„), (3) 

where f{s),g{s) and h{s) are functions of s. The La- 
grangian of the system is then written as 

71 

E([/(^"+i)-/(^")]' 

n 

+ [9{Sn+l) - giSn)? + [h{Sn+l) " KSn)f). (4) 



Then, the Euler-Lagrange equations takes the form 
"^„ [/'(s„)2 + g'{snf + h!{snf\ s„ 



= X[/(s„+i) + fisn-i) - 2/(s„)]/'(s„) 
+K[g{s„+i) +g{sn-i) - 2g(s„)]5'(s„) 
+K[h{sn+i) + h{s„-i) - 2h{sn)]h'{sn), (5) 

where A'{s) = {d/ds)A{s) with A = f,g,h. It can hence 
be seen from Eqs. (@J and ||SJ) that the geometric con- 
straint of Eqs. @, transforms the original linear lat- 
tice equations into nonlinear ones. Thus, depending on 
the choice of /, g and /i, we can obtain various nonlin- 
ear equations from the linear lattice model through the 
one-parameter geometric constraint. Equations (0) can 
be generalized to the case in which, e.g., the quantities 
Xn, yn are parametrized by two variables, (si„, S2n)- 



III. HELICAL CONSTRAINT AND THE 
SINE-LINEAR LATTICE 

We illustrate the above method by considering situa- 
tions in which the constraint functions /, g, h are of the 
form 

/(s„) = a„ cos(cs„) = a„ cos(6'„), 
gisn) = a„sin(cs„) = a„sin(6'„), 

Zn = VnSn = bnOn, (6) 

where c, am Vn are constants depending on the site index 
n and 6„ = w„/c. Inserting Eqs. © into Eq. we 
obtain the Lagrangian of the system in the form 



L^T-U = 



2)5)2 



n 

- E Y + ~ 2a„+ia„ cos(0„+i - 9n) 

n 

+ {bn+lSn+l — bnOn)'^]- 

Then, the Euler-Lagrange equations are given by 

Sn = Cn[an+ian sin(0„+i — 6n) 
-a„a„_i sin(6l„ - On-i)] 



(7) 



where 



K 



mn{al + bl)' 



(8) 



(9) 



Equations JHJ interpolate between a generalized version 
of the linear lattice 

i^n — Cn \bn+lbn()n+l + ^n^n-l^'n-l ~ '^b'^6n\ , (10) 
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and that of the sine-lattice equation |lOl 

dn — Cn[an+ia,7i Sm{9n+1 — On) 

-a„a„_i sin(6'„ - 9n-i)], (11) 

where the coefBcients of the right-hand side are all site- 
dependent. Eq. ((SJ are referred to as sine-linear-lattice 
(SLL) equations. 

IV. A PERFECT HELICAL LATTICE 

Here we restrict ourselves to the simplest situations 
in which all the coefficients are site-independent. Then, 
dropping the subscript n attached to the coefficients in 
Eq. (0 leads to 

9n = J[sin(6'„+i - 9n) - sin(6'„ - 0n-i) 

+e(0„+i+0„_i -20„)], (12) 

where 

J = Ca^, e ^ b^/a^. (13) 

The non-integrable equation (|12|l can be considered as 
relevant to helical-lattices in which the helicity is deter- 
mined by the factor e. Alternatively, Eqs. H12|l may 
be considered as a modified version of the discrete sine- 
Gordon equation '2^ in which the conventional on-site 
term sin(0„) is replaced by the first term in the right 
hand side of Eq. ifT^ . 



decrease of e, new minima appear in the potential V{5n), 
the ^th minimum appearing at 



q«(4Z-l)--^(4/-l)2-_2, ;>0. (17) 

For example, Ci « 0.2172, €2 « 0.09133, £3 « 0.05797. 
The potential V{6n) is shown in Fig. ^for these values of 
e and also for a larger magnitude, e = 0.4, when there is 
only one minimum. In the limiting case e — > 0, the SLL 
Eq. (|12|l degenerates to the SL equation with energeti- 
cally indistinguishable minima situated at Sn = 2ttI. 

As the number of minima on the potential V{6n) in- 
creases, the variety of the stable structures of the SLL 
Eq. (|12|l also increases. For e;+i < e < ei, 6n can ob- 
tain 2{l — 1) values (positive or negative) in addition to 
(5„ = 0, which is allowed for any e. 

For the particular class of equilibrium solutions ex- 
pressed by Eq. (|15l) . configuration of the nth bond, 
— 9n+i — 9n, is uot affected by the adjacent bonds, 
and one can easily create various compacton-like local- 
ized defects. For example, for e = 0.2, V{Sn) has minima 
a.t Sn — and (5„ = ±4.9063. In Fig. [3 we show sev- 
eral stable configurations of SLL at e = 0.2 and J = 1. 
Shown are (a) a kink-type structure, (b) a point defect 
in a ground state structure, (c) a zigzag periodic struc- 
ture, and (d) a point defect in the zigzag structure. The 
defects in periodic structures presented in Fig. |21can be 
regarded as discrete compactons fT7| because they are 
sharply localized and do not have exponential tails. 



A. Steady state structures 
Equation ()12|l can be rewritten in the form 

9n = J{Dn — Dn-l), 

Dn = sin{dn) + eSn, 5n = 9n+i-9n- (14) 

Equilibrium solutions to Eq. (|14|l can be found from 
Dn = Dn~i, which is particularly satisfied when 

Dn - 0, (15) 

for all n. Stable equilibria are subjected to an additional 
condition, D[^{5n) — cos((5„) + e > 0. 

To study the physical properties of Eq. ((T^ which 
represents the competition between the linearity and the 
helicity, let us consider the potential function V{5n) as- 
sociated with Dn{Sn), 

V{Sn) = l-cos{Sn) + e^, e>0. (16) 

We pay particular attention to its minimum points. For 
sufficiently large e (weak helicity), the potential V has 
only one minimum at (5„ = and the only stable state of 
the SLL is the ground state, 0„ = const. However, with 



B. Perturbation Induced Dynamical Structures 

A vibration mode can be excited on the defect in zigzag 
structure presented in Fig. |2Id). The mode is (practi- 
cally) localized on three particles and it is presented in 
Fig. 13 The existence of this mode can be easily under- 
stood. The 0th particle in Fig. |2Id) oscillates with a fre- 
quency close to the upper edge of the ground state spec- 
trum which, however, is outside of the frequency band of 
the zigzag structure. 

Due to translational invariance, the SLL model of 
Eq. H12|) admits solutions with 9n = (p for arbitrary 
if — const. Domain walls separating two domains with 
7^ V2 can also be formed. Here we consider the case 
when the difference |</3i — (^2! is not equal to the distance 
between the minima of the potential V{6n) Eq. H16|l and 
which can give rise to moving kinks. There are many pos- 
sibilities to initiate such moving kink-like structures; one 
is presented in Fig. ^for J = 1, e = 0.2. Here the kinks 
were initiated by applying initial conditions 9o(0) = 0, 
^o(O) = 3, with zero initial conditions for all other par- 
ticles. The height of the kinks is \lpi — 952! ~ 0.97r and 
it increases by increasing the initial velocity of the 0th 
particle. Here we use periodic boundary conditions to 
demonstrate that the kinks survive collisions and hence 
are relatively robust objects. 
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C. An approximate analytical solution: nonlinear 
resonant modes (NRM) 

Let us consider situation in which an atom located at 
the site makes a large excursion around one of the min- 
imum points a of the potential, while the others only 
perform small amplitude oscillatory motion. This means 
that 



|0o|»|0«| for |7^|>1. 



(18) 



As we expect a highly localized solution, we preserve the 
full nonlinearity of the equations for the central site only, 
while the other equations are linearized around the equi- 
librium position. Such a procedure leads to 

^0 = J(-2[sin(0o) + £^0] + [cos(0o) +e](ei +e-i)),(19) 



e±i = J[- [cos(6lo) + l + 2e]6l±i 
+ (l + e)e±2+sin(0o)+eeo], 



(20) 



= J(l + £)(&„+! +0„_i-20„), for|n|>2. (21) 
We now seek solutions to this set of equations in the form 

6*0 = a + Wq: with sin(Q;) + ea = 0, 

0„ = u„, for \n\ > 1. (22) 

When, in Eq. I|22|l . a — 0, all particles oscillate near 
one well of the potential V{6n), while for a 7^ the 0th 
particle oscillates in a different well. 

Inserting Eq. (EU into Eqs. pi and ^ 

and retaining solely terms linear with respect to uq and 
Unin 7^ 0), we arrive at a set of equations which are of the 
same form as those appearing in one-impurity problems 
in harmonic lattice vibrations. Thus, setting 



Vn exp(-iwi). 



(23) 



where w„ is time-independent and w is a constant, we 
reduce the above set of equations to the form 

with L^°^Vn^ J{l + e){2Vn-Vn+l-Vn-l), (24) 

with 

L'vq — J[l — cos(a)](2?;o — Vi — w_i), 
L'v±i = J[l - cos(a)](w±i - Wo), 

L'vn = for n^O,±l. (25) 

From the above equations, the ILMs we are seek- 
ing are presumed to exist within the frequency band 
w(fc)^ = 2 J(l -|- e)[l — cos(fc)] of the perfect harmonic 
lattice characterized by wavenumber k. Such a nonlinear 



localized in-band modes are referred to as (intrinsic) non- 
linear resonant modes (NRM). Then, Eqs. and 
can be handled by introducing lattice Green's functions 

g{n) = g{n,uP) 
exp{ikn) 



N ^ 2J(1 



2J(1 + e)[l - cos(fc)] - cj2 _ 



(26) 



Analytical expressions for the lattice Green's functions 
g{n) can be obtained by replacing the sum with respect 
to k by an integral over the first Brillouin zone as follows 



1 cos{nx)dx 



2 J(l + e) TT 7q y — cos(x) — ij 



2J{l + e) 



[C„(y)-z5„(y)], (27) 



where 



Cn{y) = / coa{yt)Jn{t)dt, 



Sniy) = / sm{yt)Jn{t)dt, 
Jo 

with W = 1 ; r , 

^ 2J(l + e)' 



(28) 



in which J„(t) is the Bessel function of the nth order. 
We note that the quantity 



C'oiy) = 



1 



with So{y) for < y < 1, 



(29) 



represents the density of states of the band. In terms of 
the gnS so obtained, Eqs. H24|l and (|25|l are rewritten as 

Un = J [1 - cos(fc)] [g{n){2vo - ui - w_i) 
+g{n - l){vi ~ vo) + g{n + l)(w_i - vq)] . (30) 

We pay particular attention to an s-like mode having the 
symmetry property 



and use the identity relation. 



(31) 



J{l + e)[2g{n)-g{n+l)~g{n~l)] 

-uj^g{n) = A„, (32) 



to obtain 



where 



= A[w^g(n) + A„] [vq - vi), 



1 — cos(a) 
1 + e ' 



(33) 



(34) 
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and A„ are Kronecker's delta. 

In terms of the dimensionless frequency 
^ = J(l + e), an equation giving the eigenfre- 
quency of the NRM is obtained from Eqs. lt?7|l . (|^ . 
and H33|l as follows: 

1 - A - AC - iXeCoiy) = 0. (35) 

The dimensionless squared eigenfrequency ^ is therefore 
obtained as 

The localization of the NRM is given by the equation 

Vn = Xuj'^g{n){vQ - vi), n ^ 0. (37) 

As mentioned before, here the localized mode around the 
local minimum point of the potential function is of reso- 
nant type, the eigenfrequency which appears close to bot- 
tom of the frequency band uj{k)'^ = 2J(1 -I- e)[l — cos(A;)]. 
Apart from the factor a, which is the position of the lo- 
cal minimum, the localization properties of the NRM are 
essentially different from those of the ILM in that it ex- 
hibits oscillatory slow decay, contrary to what is the case 
for the ILM. This point will be examined in detail below. 

We have attempted to excite NRMs in SLL of Eqs. 
(|12|l . Such a numerical calculation is presented in Fig. |S1 
Here we set J = 1, e = 0.2 and for the initial conditions, 
^0(0) = 4.88, 0o{O) — 0, with zero initial conditions for 
all other particles. For this choice of e, the potential 
V{Sn) has minima at (5„ — and (5„ « ±4.9063. One 
can see that the 0th particle oscillates near the minimum 
of V{6n) situated at Sn ~ 4.9063, while the remaining 
particles oscillate near the minimum at 5„ = 0. Due to 
the rather strong interaction with the linear excitations 
of the lattice, the energy of the NRM is rapidly imparted 
to the lattice (after a few oscillations). The lifetime of 
the excitation does not increase much, either by changing 
the initial deviation of the 0th particle from the potential 
minimum or by changing e. 

V. INCLUSION OF THE ON-SITE POTENTIAL 

A. Derivation of the model 

It is of interest to see what happens if we include on- 
site potentials in our original lattice. This is done by 
generalizing Eq. ^ in the form 

n 

n 

+ {zn+i-z^f + ^xl+'^yl], (38) 



where ki and K2 are constants. Applying the helical con- 
straint, Eq. ((HJ leads to 

n 

- ^ y [a-n+i + On ~ 2a„+ia„ cos(6'„+i - On) 

n 

+ {bn+lOn+l — hnOnY 

+ ^^^^^^^cos[l-cos(20„)]]- (39) 

It is seen that such a procedure adds a sine-Gordon type 
on-site potential 1 — cos(26'„) to the Lagrangian. The 
equations of motion for 0„ then assume the form 

6n — Ln [an+lO-n sin(0„+i — On) 

-a„a„_i sin(6'„ - On-i) 
-Xalsm{20n)] 

-^Ln\bn+lbnOn+l + ^n&n-l^n-l 

-{hn+lhn + bnhn-l)On], (40) 

where A = {n2 — ki)/2. When all the coefficients are site- 
independent, the above equations reduce to 

On = J[sin(6'„+i - On) - sin(6'„ - On-i) 

-X sin(20„) -I- e{On+i + On-i - 20„)] . (41) 

The presence of the on-site potential here can be spec- 
ulated to induce the existence of localized modes with 
longer lifetime than that of the case of Eq. (|12(l . In the 
following, we demonstrate the existence of ILMs in the 
model of Eq. |(1T|| . 

B. Intrinsic localized modes 

The linear spectrum of the SLL with the on-site po- 
tential, Eq. (|41|l . is given by 

a;2 = J[4(H-e)sin2(fc/2) -f 2A], (42) 

with 

max — v/J[4(l + e) + 2A]. (43) 

When Wmax/2 < Ci^min, 1-6., A > (2/3)(l + e), one can 
look for an ILM oscillating with frequency tj,nax/2 < < 
^min so that all higher harmonics Zw, with an integer 
I > 1, lie above Wmax- In this situation the ILM would 
not interact with the lattice preserving its identity. 

1. ILM with all particles in one potential well 

As an example, we take J = 1, e = 0.2, and A = 5, 
with Wmin = \/T0 ~ 3.162, and Wmax = V14.8 « 3.847 
so that the necessary condition of the existence of ILM, 
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'^max < Zwmin, IS fulfilled. A large- amplitude localized 
mode can be excited by choosing the initial deviation 
and/or initial velocity of a particle. As an example, we 
initialize the 0th particle with 6*0(0) = 1.5, Oo{0) = 
with zero initial conditions for other particles and, after 
some stabilization period, a steady oscillatory motion is 
observed (see Fig. (HI . One can see that the ILM is highly 
localized. The mode has a frequency lo « 2.15, which is 
considerably lower than the bottom edge of the linear 
spectrum, Wmin, with all higher harmonics lying above 
the upper edge, Wmax- 

The ILM can be expressed analytically assuming that 

9o{t) Ai sm{ujt) + A3 sin(3wt), 
6»i(t) = e-i{t) = Bi sin(wt), 

0„(t) = for \n\ > 1, (44) 

with A3 <C Ai and Bi <C Ai. Parameters of the approx- 
imate solution Eq. (|44|l can be expressed in terms of the 
ILM amphtude, Ai: 



J 



2{1 + e + X) - hi + 4X)Al 



A. 



(1 + 4A)A? 



12 [9^ - 2(1 + e A)] ' 
(1 + e)A^ - iA? 



2(1 



A)-^ 



(45) 



Let us analyze the obtained solution. First of all we 
note that the ILM frequency, uj, cannot be greater than 
'^min, i-e., we can only look for a mode with frequency 
below ojmin [see Eq. H43(l ]. which is possible only for ILA4s 
with sufficiently large amplitude. 



Ai > 



8(1 



1 + 4A 



(46) 



For example, for e — 0.2, and A = 5, one has Ai > 0.68. 
We have confirmed numerically that, for this choice of 
parameters, our solution can be used for Ai ~ 1. 

For larger amplitudes the ILM is highly localized (see 
Fig. O but the solution Eq. igJJ), becomes inac- 
curate because it takes into account only cubic anhar- 
monicity. The ILM in Fig. Qhas amplitude Ai « 3.71 
and frequency lo ~ 0.607, which is significantly lower the 
bottom edge of the linear spectrum, Wmin ~ 3.162. The 
mode radiates energy extremely slowly due to the interac- 
tion of higher harmonics with the linear spectrum. Note 
the difference in the ordinate scale for the middle panel. 

For small ILM amplitudes, approaching the allowed 
limit, Eq. (|46|) . the solution profile becomes wider and 
our assumption that it is localized on three particles be- 
comes invalid. Estimation for the limiting amplitude, 
Eq. 1)46(1 . is also valid only for a highly localized ILM. 
In fact, it is possible to excite an ILM with a very small 
amplitude but the width of such ILM increases consider- 
ably and its frequency approaches LOmin from below. The 



smooth transformation of the ILM from a highly localized 
(as the amplitude decreases) into a very broad one is not 
surprising because for smooth solutions, (0„+i — 0„) <C 1, 
one has sin(^„4_i — 9n) ~ (^n+i — ^n) and the SLL with 
the on-site potential, Eq. H41() . can be approximated by 
the well-known continuum sine-Gordon equation. The 
ILM is, then, nothing but the analog of breather solution 
to the sine-Gordon equation. Using the Lorentz invari- 
ance of the sine-Gordon equation one can create a moving 
small-amplitude breather. 

We also note that, when e ^ 1 and A ^ 1, the first two 
terms in the right hand side of Eq. H41|) can be neglected. 
In this situation, when e and A are of the same order of 
magnitude, our solution, Eqs. (|44|l and H45(l . describes 
a highly localized ILM in the Frenkel-Kontorova model 

iai. 



2. ILM with one particle in a different potential well 

The SLL lattice with on-site potential supports local- 
ized modes of another type, when one particle is trapped 
in a potential well different from the well occupied by 
the other particles. To demonstrate this mode, we take 
J = 1, e = 0.05, and A = 1, with Wmin = V2 « 1.414, 
and Winax — %/ 6.2 w 2.490 so that the necessary con- 
dition, Wmax < 2a;niin, IS fulfilled. Setting 6*0(0) = 2.9, 
^o(O) = with zero initial conditions for other particles, 
after some time a steady oscillatory motion presented 
in Fig. IHlwas observed. The 0th particle was found to 
oscillate with to « 1.35 which is below the linear vibra- 
tion band but all higher harmonics are above the band. 
One can see that the ILM is not highly localized. The 
0th particle oscillates with the amplitude A{9o) ~ 0.6 
near the coordinate shifted by roughly 2n with respect 
to the neutral position of all other particles [shown in 
Fig. Elb)]. Its nearest neighbors oscillate with nearly the 
same amplitude, A(6±i) « 0.52, and the amplitudes de- 
crease rather slowly with deviation from the 0th particle 
[see Fig. IHIa)]. 



VI. 



CONCLUDING REMARKS 



In this paper, we have outlined a general scheme to 
apply geometrical constraints to a simplified version of 
the 3D harmonic lattice with or without an on-site po- 
tential to obtain nonlinear dynamical lattice equations. 
Depending on the type of the constraints, we arrive at 
various nonlinear equations, being interested in the co- 
herent structures that arise in them. 

We then studied, more specifically, a helical constraint 
by which the original 3D harmonic lattice was shown to 
reduce to ID sine-plus-linear lattice equations or helical 
lattice equations. This formulation revealed several in- 
teresting features: 

1. It provides us with a systematic derivation of the 
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sine-latticc equations which were derived heuristi- 
cally in earlier works. 

2. Physically, applying the helical constraint amounts 
to transforming the original harmonic lattices (pos- 
sessing the hard potential) to a nonlinear dynam- 
ical lattice that supports resonant and multi-kink 
modes. Seeking other kinds of nonlinear modes will 
be an interesting topic for future study. 

3. The method developed here may be relevant to the 
study of nonlinear excitations in biomolecules with 
helical structure such as DNA and proteins. There 
are several areas of condensed matter physics where 
the existence of a complex energy landscape with 
a number of local minima in the potential function 
is presumed to play a crucial role in determining 



their properties. Such examples can, for instance, 
be found in the study of glasses and proteins. Ap- 
plication of the techniques employed here to these 
problems would be another topic warranting fur- 
ther investigation. 

4. Finally, another relevant generalization would be 
to consider geometric constraints different than the 
helical one, and derive and study the ensuing (re- 
duced) dynamical models. 

These topics are currently under study and will be re- 
ported elsewhere. 
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FIG. 1: The potential V{&n) of Eq. GgJ is shown as a function 
of 5n for different magnitudes of e. Curves 1 to 4 correspond 
to e = 0.4, ei « 0.2172, ea « 0.09133, and eg ~ 0.05797. 
The magnitudes e;, at which /th minimum of the potential 
function appears, are given by Eq. H17^ . 




FIG. 2: Examples of stable equilibrium structures in SLL with 
J = 1, £ = 0.2. (a) Kink-like structure, (b) point defect in the 
ground state structure, (c) zigzag structure, (d) point defect 
in the zigzag structure. 
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FIG. 3: Vibration mode excited for zigzag metastable struc- 
ture with a defect present in Fig. Eld). The mode is localized 
on (practically only) three particles. J — 1, e — 0.2. 
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FIG. 4: Moving kinks initiated by the initial condition 6o{0) — 
0, 6o{0) ~ 3, and zero initial conditions for all other particles, 
presented by the set of plots 9 vs. n for different times. The 
height of the kinks is \ipi — (p2\ ~ O.Qtt and it increases with 
increase in initial velocity of 0th particle. Periodic boundary 
conditions are used here to demonstrate the kinks' robustness 
against collisions. J — 1, e = 0.2. 
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FIG. 5: An attempt to excite a NRM in a SLL lattice Eqs. 
l(T^ by means of initial conditions ^o(O) = 4.88, ^o(O) = 0, 
with zero initial conditions for all other particles. J = 1, e = 
0.2. The 0th particle oscillates near the minimum of V{S„) 
situated at 5„ ~ 4.9063, while other particles oscillate near 
the minimum at 5„ = 0. Due to the rather strong interaction 
with the lattice, the energy of excitation is "dissipated" to 
the lattice after only a few oscillations. 
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FIG. 6: ILM with all particles oscillating in one well of the 
on-site potential. The mode has amplitude Ai ~ 1.37 and 
frequency ui ~ 2.15, which is considerably lower than the 
bottom edge of the linear spectrum, Wmin « 3.162 with all 
higher harmonics lying above the upper edge, 3.847. 
J = 1, £ = 0.2, and A = 5. 
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FIG. 7: Same as in Fig. |S|but for larger amplitude, Ai ~ 
3.71. The ILM has frequency lj ~ 0.607 which is significantly 
lower than the bottom edge of the hnear spectrum, LOmin ~ 
3.162. The mode radiates energy extremely slowly due to 
the interaction of higher harmonics with the linear spectrum. 
Note the difference in the ordinate scale for the middle panel. 
J = 1, e = 0.2, and A = 5. 
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FIG. 8: ILM with one particle oscillating in a potential well 

different from that where other particles are located, (a) Plots 
9 vs. n for different times showing the dynamics of parti- 
cles (the 0th particle is not shown hero because its coordi- 
nate differs from coordinates of other particles significantly, 
by roughly 2tt). (b) Dynamics of 0th particle. The ILM is 
rather broad. The 0th particle oscillates with the amplitude 
A{do) ~ 0.6 while its nearest neighbors have nearly the same 
amplitude, A{9±i) « 0.52, and the amplitudes decrease rather 
slowly with deviation from 0th particle. J = 1, e = 0.05, and 
A = 1. 



